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' 

' The exhaustive classification of stationary incompressible flows with constant total 

' pressure of ideal infinitely electrically conducting fluid is given. By introduction 

of curvilinear coordinates based on streamlines and magnetic lines of the flow the 
■ system of magnetohydrodynamics (MHD) equations is reduced to a nonlinear vector 

wave equation extended by the incompressibility condition in a form of a generalized 

1 . Cauchy integral. For flows with constant total pressure the wave equation is explicitly 
j—t^ . integrated, whereas the incompressibility condition is reduced to a scalar equation 

^ . for functions, depending on different sets of variables. The central difficulty of the 

i investigation is the separation of variables in the scalar equation, and integration of 

the resulting overdetermined systems of nonlinear partially differential equations. The 
I canonical representatives of all possible types of solutions together with equivalence 



Ph. 



> 



transformations, that extend the canonical set to the whole amount of solutions are 
found. 



Introduction 

The mathematical model of ideal magnetohydrodynamics (MHD) describes macroscopic 
motions of infinitely electrically conducting fluid under the action of internal pressure, 
magnetic and inertia forces. The area of applicability of the model is extremely large: 
from problems of magnetic confinement to astrophysics. The popularity of the model is 
explained by its relative simplicity, by the wealth of its mathematical content, and by the 
variety of described physical effects. In most applications it is required to have solutions 
that model essentially three-dimensional motions of fluid. At that, the numerical or 
analytical analysis is complicated by the nonlinearity of equations. In this respect exact 
solutions that describe three-dimensional flows of conducting fluid are of interest. 

Advancement in the construction of exact solutions can be reached by introduction 
of special curvilinear system of coordinates [4], determined by the geometry of 
streamlines and magnetic lines. In this form MHD equations are partially integrated 
and reduced to a vector wave equation extended by the incompressibility condition in 
the from of the generalized Cauchy integral. Any solution to the equations immediately 
gives explicit description of magnetic lines and streamlines of the flow. 
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This work is devoted to complete description of ideal electrically conducting fluid 
flows with constant total pressure p+ ||Bp = const. For this class of solutions the 
vector wave equation is integrated explicitly. Substitution of the obtained representation 
of solution into the Cauchy integral yields a nonlinear scalar equation, that contains 
functions, depending on different sets of arguments. Separation of variables in the scalar 
equation, and subsequent integration of obtained overdetermined systems of partially 
differential equation is the main difficulty of the research. The significant simplification 
of calculations is reached by the systematic use of equivalence transformations, which 
allow one to discard insufficient arbitrary functions and constants. The separation 
of variables performed in Theorem [1] produces a set of 9 overdetermined systems of 
equations. Integration of each system with the use of equivalence transformations leads 
to the formulated in Theorem [2] main result of the paper: the complete list of canonical 
forms of solutions, describing MHD fiows with constant total pressure. It is shown that 
in all such fiows Maxwellian surfaces bearing streamlines and magnetic lines of the fiow 
belong to a class of translational surfaces, i.e. are swept by the parallel translation of 
some 3D curve, called generatrix, along another 3D curve called directrix. In Theorem [3] 
it is proved that all solutions belong to either of two main classes: fiows with cylindrical 
Maxwellian surfaces, or fiows with Maxwellian surfaces having a common directrix. 

Class of stationary ideal MHD fiows with constant total pressure was earlier 
investigated in J4]. Here the curvilinear system of coordinates was introduced, general 
properties of solutions were investigated, some exact solutions describing motions with 
constant total pressure were obtained. In [6] with the use of a similar curvilinear 
coordinates there were obtained exact solutions, where pressure is constant not in the 
whole area occupied by the fiow, but only on toroidal Maxwellian surfaces. Classes 
of solutions with constant total pressure in case of non-stationary fiows were obtained 
in [5]. In particular, solutions describing fiows with knotted or toroidal Maxwellian 
surfaces were constructed. Note, that problem of complete description of such solutions 
in non-stationary case is technically more complicated and is not solved yet. 

1. Preliminary information 

Equations of ideal magnetohydrodynamics (MHD), describing stationary fiows of ideal 
infinitely conducting fiuid with frozen-in magnetic field have the following dimensionless 
form: 

p(u ■ V)u + B X rot B + Vp = 0, rot (u x B) = 0, 

(1.1) 

divu = 0, divB = 0. 

Here u is the velocity, B is the magnetic field, p is the density, p is the pressure. Density 
p is assumed to conserve along streamlines: u ■ Vp = 0. The following functions are 
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introduced: 

v = uv^, P = p+^B^ (1.2) 

Function P will be referred to as the total pressure. With these functions system (11.11) 
can be rewritten in the following form: 

(v V)v-(B-V)B + VP = 0, 

(B- V)v = (v V)B, (1.3) 
divv = 0, divB = 0. 

Equation for density p splits off from the system ( 11. 3p and can be solved separately. 
Equations (11.31) form an overdetermined compatible system of 8 equations for 7 sought 
functions. 

Except for the obvious symmetry group of transformations consisting of shifts along 
the coordinate axes, rotations, dilatations and shift of the complete pressure, system 
of equations (11.31) admits the following transformation of velocity and magnetic vector 
fields [HE]: 

V = B sinh /((/?) + V cosh /((/?), B = B cosh /((/?)+ v sinh /(</?). (1.4) 

Here cp is an arbitrary function that conserves along streamlines and magnetic lines of 
the flow: 

(vV)<^ = 0, (B-V)^ = 0. (1.5) 

The symmetry group of system fll.3p will be used for simplification of the form of 
obtained solutions. 

2. Curvilinear system of coordinates 

System of equations (I1.3P has more symmetric form in terms of the following vector 
fields, also known as Elsaasser variables [2]: 

a = v-B, b = v + B (2.1) 

or 

v = l(b + a), B = i(b-a). (2.2) 
In terms of these vector fields system (II. 3p is transformed to 

(a-V)b + VP = 0, (2.3a) 

(a- V)b = (b- V)a, (2.35) 

diva = 0, divb = 0. (2.3c) 
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Note that equation (12.36P is equivalent to the equahty to zero of the commutator 
of vector fields a and b: 

[a,b] := (a- V)b- (b- V)a = 0. (2.4) 

Under the assumption of linear independence of these vector fields, the commutativity 
condition (12 ■4p implies that the fields can be chosen as basic ones for a certain curvilinear 
system of coordinates [7]. Indeed, let us observe a curvilinear system of coordinates 
[k^, A;^, k^) such that k^- and fc^-coordinate lines coincide with integral curves of vector 
fields a and b accordingly. In other words, in the space the following non-degenerate 
mapping is introduced 

x = x(k), det(^^^^0 (2.5) 

such that 

9x , 9x , . 

In this representation equation (12. 4p is satisfied automatically as a compatibility 
condition of equations (12. 6p for unknown vector x(k). The remaining equations (I2.3ap . 
(12.3 cp lead to the following system of equations for unknown functions x = x(k), 
P = P(k) (details can be found in [4]): 

9^x dP f dx dx\ dP / (9x &x\ dP ( (9x 9x 

)=U, 

(2.7) 



dk^dk'^ ' dk^ V dk'^ ^ dk^J ^ dk^ \ dk^ ^ dk^J ^ dk^ \ dk^ ^ dk^ J °' 



5x 



dP 

^ dk^ 


r dx 


\dk^ 


dk^ ^ 


dx\ 


dk^j 



1. 



dk^ 

Henceforth " ■ ", and " x " denote the Euclidean dot and vector products in M"^ 
correspondingly. The first equation in (12. 7p can be treated as a nonlinear vector 
wave equation for components of vector x. The second equation in (12.71) represents 
an incompressibility condition of both fluid an magnetic field. System of equations (12. 7p 
is equivalent to the system ( I2.3ap -( j2.3cp on solutions with non-collinear vector fields a 
and b. 

Parametrization of the solution in terms of the dependence x(k) gives the explicit 
description of the geometry of the fluid motion. For example, important in applications 
Maxwellian surfaces bearing magnetic field lines and streamlines of the flow, in this 
parametrization are given by x = x(/c^, A;^, /cq) with arbitrary parameters k^, fc^, and 
fixed parameter k^ = /cq. Each Maxwellian surface can be taken as an infinitely 
conducting wall bounding the area of the flow. 

Note that system of equations (12. 7p is invariant under the following change of 
parameters 

~k' = ip{k')k\ P = ^, P = P (2.8) 

(flk-^j 
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with an arbitrary function y?. This symmetry follows from the transformation fll.4p . 
admitted by the initial system (11. 3p . Besides, system of equations (12. 7p is invariant 
under transformations 

fci = A;^ + ~k^ = e + x{k^) (2.9) 

with arbitrary functions and x- what follows formulae (12. 8p . (12. 9p will be used as 
equivalence transformations on the set of solutions of equations (12. 7p . 



3. Flows with constant total pressure 

System (12. 7p possesses a class of solutions where the total pressure P is constant over 
the whole area occupied by the flow: 

P = const. 

In this case the first equation of (12. 7p is reduced to a linear one, and can be integrated 
explicitly as 

^^ = ^ ^ = a{k\e) + r{k\e). (3.1) 

Here cr and r are vector functions depending on different independent variables. 
From (13. ip it immediately follows that Maxwellian surfaces on this class of solutions 
are translational surfaces. They are spanned by parallel translation of the curve 
X = (T{k^, /cq) along the curve x = t(/c^, /cq) with a fixed parameter k^ = /cq. 

Substitution of representations (13.11) into the second of equations (12.71) gives the 
restriction for vectors cr and r: 



da 



dr ( dcr dr 



dk"^ \ dk^ dk^ 



1. (3.2) 



Now the description of flows with constant total pressure is reduced to a separation of 
variables k^ and k"^ in equation (13.21) and integration of obtained overdetermined systems 
of equations. 

Henceforth the following notations are adopted. The standard Cartesian basis of 
orthonormal constant vectors in R'^ is denoted by e^, e^, e'^. The lower index i denotes 
the derivative with respect to k^ (e.g., cri = dcr/dk^). Bold Roman or Greek letters 
denote vectors. 



4. Separation of variables 

Note that equations (13. 2p can be written as the equality to unity of the dot product of 
two vectors in R^, such that each vector depends on its own set of variables: 

A-B = l, (4.1) 
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where 

A = (<T3 X CTi, (Tif, B = (t2, T2 X T^f (4.2) 

Here A = A{k^^k^)^ B = B(A;^,A;^), symbol "T" denotes the transpose of a matrix 
or a vector. Separation of variables and k"^ in equation (14. ip is managed by the 
following Lemma, valid for general n-dimensional vectors A(x), B(?/) depending on 
different variables x and y 

Lemma 1. Equation i4.1\) is satisfied for vectors 

A = (ai(x), . . . , a"(x))^, B = (^^(y), . . . , }f{y)f 

identically on variables x and y if and only if there exists a number s {0 < s < n — 1), 
constant n x s matrix C of rank s, constant vector d G M" and a set of s non-constant 
linearly independent functions , forming a vector u{x), such that, the following equations 
are satisfied: 

A = Cu{x) + d; B^C = 0, B^d = 1. (4.3) 

Proof. Let the linear span of vectors A(x) for various admissible values of x has a 
dimension k. This implies that there exist exactly k various values x^, i = 1, . . . ,k, for 
which constant vectors A* = A(x*) are linearly independent. An arbitrary vector A(x) 
can be represented as a linear combination of the constant vectors as 

k 

A{x) = J2u'{x)A\ (4.4) 

1=1 

By virtue of (14. 4p one has u^{x^) = 5^^ with the Kronecker delta in the right-hand side. 
This implies linear independence of functions m*(x). For the compatibility of equations 
(14. ip one function among should be constant. Suppose, that this function has index 
A;, i. e. = 1. By arranging vectors A^,...,A'^~^ in a matrix C with n columns 
and s = k — 1 rows and denoting d = A'^ one obtains a representation of vector A(x), 
specified in the formulation of the Lemma. Equations for vector B(y) follows from the 
substitution of the representation of vector A into equation (14. ip and splitting with 
respect to linearly independent functions u^{x). □ 

In what follows Lemma [1] is used as the main tool for separation of variables in 
equation (14.11) with vectors A, and B determined by (14. 2p . At that, parameters k^ and 
k"^ play the role of variables x, and y. Matrix C and vector d are supposed to depend 
arbitrarily on variable k^. Only cases s = 0, 1,2 will be successively investigated. By 
virtue of the symmetry of equation (14.10 with respect to the change of variables k^ <^ k"^, 
cr ^ T, cases s = 3, 4, 5 are equivalent to the observed ones. 

In what follows, the following equivalence transformations of equation (14.11) will be 
used for simplification of the resulting equations and solutions: 



• a = (T — e{k^) , T = T + e{k^) with arbitrary vector e; 

• Transformations (12. 8p and (12.91) . 

The analysis is preceded by the description of two characteristic cases of integration 
of equation ( 13. 2p . 



5. Case of cylindrical Maxwellian surfaces 

Suppose that vector cr depends linearly on variable k^: 
a- = k^r){k^). 

Substitution of this representation into equation (13.21) gives 

fcVs ■ (77(fc3) X V(A;3)) + r7(P) ■ (t2 x rg) = 1 

Hereafter prime denotes the derivative with respect to the single variable k^. Splitting 
on k^ produces a system of equations for vector r: 



T2 • {ri'in X r){k')) = 0, r){k') ■ {r^ x rg) = 1. 



(5.1) 



The first of equations (15.11) is satisfied identically if 77 = const. In this case without loss 
of generality one can assume rj = e^. The exact solution can be written as 



x = A;^ei + r(r,r 



^1 ^1 
^1 ^1 



1. 



(5.2) 



For a non-constant vector 77 the first of equations (15.11) gives 

r = u{k\ e)r}'{e) + v{k\ k'')'n{e) + di^). 

By substituting this representation into (15. ip after obvious transformations one obtains 
an equation for function u in the form 

du 



dk'' 



Integration with respect to k'^ leads to the following exact solution: 

X = k^'qi.k^) + u{k^, k^)'n\k^) + t;(A;2, fc3)^(P) + ^(^3^^ (5 3) 

Here rjlk^), and d{k^) are arbitrary vectors, f (A;^, k^) is an arbitrary function. Function 
u{k'^, k^) is determined by the solution of the quadratic equation 

''■^'^'"^"^^■^ ' {ri-{rj'.5'))u = k\ 



u 



(5.4) 



Due to the linear dependence of vector x on parameter fc^, Maxwellian surfaces in this 
solution are exhausted by cylinders with the generatrix directed along vector /3. At 
that, the directrix of the cylinder is defined with a significant functional arbitrariness 
by formulae above. 
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6. Maxwellian surfaces with a fixed directing curve 



Equation (13. 2p also possesses the foUowing class of solutions: 
^ = (T{X{k\e)) + T{k\k^), 



(6.1a) 



a{\{k\k^))-a{k^) = k\ 



(6.16) 



r2 X Ts = cx{k^). 



(6.1c) 



Here cr(A) and cx{k^) are arbitrary vectors. Scalar equation ( ]6.16[) implicitly defines 
function \{k^ , k^). Equations (16. Id) with fixed right-hand side form a closed system of 
differential equations for three components of vector r. 

Remark 1. Solution (I6.1ap - (l6.1cp can be generalized in the following way: If one of 
components of vector cr is identically zero, then the corresponding component of vector 
ex can be assumed to be a function of variables /c^, and k^ defined by vector r from the 
same component of equation (I6.1cp . In this case the total number of equations (16. Id) 
for vector r reduces. 

In what follows we will refer to solution (I6.1al) -( l6.1d) taking into account Remark [H 
In this class of solutions Maxwellian surfaces have common generatrix given by equation 
X = cr(A). Parametrization of this curve varies from one Maxwellian surface to another, 
however the generatrix itself remains unchanged. 

7. Separation of variables (13. 2 p 

By using Lemma [H, after partial integration and simplification of the equations, one 
obtains the following 

Theorem 1. Equation Ii3.2\) after the separation of variables can be brought by the 
equivalence transformation to one of the following non- equivalent forms: 



s = 








(7.1) 



s = 



1 




CTi = U{k^, k^)e + k; 
0, K ■ (t2 X Ts) = 1. 



(7.2) 



CTg X cTi = 1/ X e, CTi = U(k^, k^)e; 
. (e X T2) = 1, T3 ■ (e X T2) = 0. 



(7.3) 



fx ■ {n X T2) = 0, T3 ■ (k X T2) = 1. 



0-3X0-1= 0, (Ti = U\k\ + f/2(^l^ f.3y ^ ^. 



[T2 X T3) 



0, 



(T2 X r3 



K ■ (T2 X T3) 



CTq X CTi 



el X e2, 0-1 = U\k\ k^)e^ + U\k\ k^)e^- 



T2 ■ (el X e^) 
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(7.4) 

(7.5) 
(7.6) 



1, el ■ (t2 X T3) =0, e'- (t2 x T3) = 0. 
cr3 X cTi = (/"^{k^, k^)e x k, (Ti = U^{k^, k^)e + n, 
T2- {ex k) = 0, e ■ (t2 x T3) = 0, k • (ts x T3) = 1. 

0-3X0-1 = {U^{k\ k'^)n + «/) X e, CTi = [/i(A;i, A;3)e; 
r2-(/xxe) = 0, T2-(i^xe) = l, e ■ (t2 x T3) = 0. 
0-3x0-1 = {U\k\ k^)fi^ + U\k\ k^)^i^) X 0-1 = k; 

T2 ■ (/xl X k) = 0, T2 • (/X^ X k) = 0, K ■ (t2 X Ts) = 1. 

In these equations the sought vectors are cr{k^,k^), and T{k'^,k^); their lower indices 
i denote derivatives with respect to k\ All the remaining vectors depend arbitrarily on 
variable k^, and are linearly independent for every value of k^ . Arbitrary functions 
U , f/i, U"^ , and 1 are linearly independent as functions on variable k^ over the whole 
domain of each system. 



(7.7) 



(7.8) 



(7.9) 



Proof of Theorem [T] can be found in Appendix A Each of the systems f l7.ip -( 17^ 
is a nonhnear overdetermined system of differential equations with sought functions cr, 
and T. Each system possesses its own group of equivalence transformations. The group 
will be utilized to discard the inessential arbitrariness of arbitrary functions and vectors 
enumerated in the formulation of Theorem [H 

In the following two sections the investigation of compatibility conditions and 
integration of each of obtained in Theorem [1] systems of equations is performed. The 
equivalence transformations of each system are used in order to bring its solution to a 
certain canonical form. The result of computations is formulated in Theorem [2l The 
general standard form of the solution is given in Theorem [3l 

8. Integration of equations for s = and s = 1 

The analysis is started from system (17. ip . Integration on k^ of the second equation of 
( 17. ip gives 

o- = k^u{e) + 5{k^). 
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Substitution of this result into the first equation of 07.11] leads to 

As before, prime denotes the derivative with respect to the only argument k^. Splitting 
with respect to the independent variable k^ and integration of the subsequent equations 
produce 

ly = a{k^)r], 6 = l3{k^)v, rj = const. 

with an arbitrary function a{k^) and a constant vector rj. By the action of 
transformations (12.81) and (12. 9p it is always possible to fix constants as a = 1, /3 = 0. 
Thus, vector cr has the form 

cr = k^r). (8.1) 

The last of equations (17. ip after the substitution of a from (18. ip can be written as 
follows: 

77 ■ (r2 X T3) = 1. (8.2) 

The obtained solution coincides with (15. 2p . 

For s = 1 equations (17. 2p are considered. From the first of equations (17.20 it follows 
that vectors cr^ and cri are proportional which, in turn, imply 

d\ 

^ = $(A), \ = \ik\e), —^0. 

Integration of the second of equations ( I7.2p on k^ reduces the problem of determination 
of the vector field cr(/c^, k^) to the solution of the following functional equation 

^iX{k\ k^)) = u{k\ k^)e{k^) + k^K,{e) + d{k^). (8.3) 

Here du/dk^ = U . It is required to find all possible functions A(A;^, /c^), and vector fields 
$(A), such that equation (18.30 is satisfied identically on variables k^ and k^ for some 
functions u{k^,k^), d'^u/dk^'^ 7^ and vector fields e{k^), i^{k^), <5(P). 
Differentiation of the equality (18. 3p with respect to k^ yields 

$'(A)Ai = uie + K. 

Here, as usual, the lower index denotes the derivative with respect to k^. Vector 
multiplication of this equations on e, and subsequent dot multiplication on k gives 
the following classifying equation: 

Ai($' X e) ■ K = 0. 

By definition, Ai 7^ 0, hence the mixed product is equal to zero: 

■ (e X k) = 0. 

Note, that vectors and e x k depend on different sets of variables, which allows 
application of a modification of Lemma [T] for the separation of variables in this equation. 
Due to the linear independence of vectors e and k, the following cases are possible: 
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((a)) $(A) = /(A)ei, at that ■ (e x k) = 0; 

((b)) $(A) = /(A)ei + ^(A)e2, at that ■ (e x k) = 0, ■ (e x k) = 0. 
Here and are orthonormal constant vectors. 



In case (a) one has 

Substitution of obtained representations into equation (18. 3p yields 

f(\)e^ = {u{k\ P) + k^w{k'^))e{k'^) + k\{e)e^ + (5(P). 

Vector muhiphcation of this equation to and sphtting on k^ gives e x = 0, 
(5 X = 0. Hence, vectors e and 5 are also proportional to vector e^. Substitution of 
vectors n and e into the second of equations (17.31) leads to the contradiction. 
Now, case (b) is observed. Here 

K = v{k^)e^ + w{k^)e'^, e = m{k^)e^ + n{k^)e'^. 

Form equation (18.31) it follows that vector 6 is in fact a linear combination of vectors 
and e^. Thus, the equation follows: 

cT = fiX)e' + giX)e\ (8.4) 

Function X{k^, k^) satisfies an implicit equation obtained by the elimination of function 
u{k^, k^) from the following two relations: 

/(A) = u(k\ k^)m(k^) + k^v(k^) + a(k^), 

(8.5) 

^(A) = u{k\ k^)n{k^) + k^w{k^) + (3{k'^). 

Here m, n, v, w, a, P are arbitrary functions of k^ such that 

A = nv — mw ^ 0. (8-6) 

Vector cr determined by equation (18.41) . (18.51) has the form 

(Ti = u{k\ k^){m{k^)e^ + ri(P)e2) + v{k^)e^ + w{k^)e^. 

Substitution into equation cri ■ {t2 x T3) = 1 and splitting with respect to k^ lead to 
the equation 

(me^ + ne^) ■ (t2 x T3) = 0, 

(8.7) 

(f + we^) ■ (t2 X T3) = 1. 

Vector T has the following components in the orthonormal basis of vectors e^, and 
e^: 
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Equations f l8.7p can be written as follows: 



m 


n 







V 


w 





rl 


rl 


'2 


= 0, 


rl 


rl 


'2 
















rl 


rl 


^3 




rl 


rl 


^3 



(8.8) 



As before, lower indices denote derivatives on variables fc^ and k^. By virtue of inequality 
f l8.6p . equations fl8.8p can be represented in the form of two linear equations for functions 
r^(/c^, k^) and T'^{k'^, k^) provided function T^{k'^^ k^) is arbitrarily fixed: 



rl 


'2 


m 


rl 


'2 


n 


rl 


^3 


" A' 


rl 


^3 


^ A 



nv — mw. 



(8.9) 



(8.10) 

is determined by the implicit 



A = 

i ri A' 

Finally, the solution takes the form 

X = (/(A) + T')e' + ((7(A) + T^)e^ + T^e\ 

where /(A) and g{\) are arbitrary functions; X{k^,k^) 
equation that is obtained by the elimination of function u{k^,k^) in relations fl8.5p : 
T^{k'^,k^) and T'^{k'^,k^) are solutions of linear equations (18.90 with arbitrarily fixed 
function T^{k'^, k^) and arbitrary functions m{k^), n{k^), v{k^), w{k^), taken form (18. 5p . 
This solution is simplified below by the action of equivalence transformations. 

Let m = 0. By using the equivalence transformations (12.80 and (12.90 one can make 
f = 1, a = in the first of equations (18. 5p . According to the obtained equation, by 
virtue of the arbitrariness of function /, one can assume A = k^. At that, the second of 
equations (18.50 does not lead to additional restrictions to functions g and A due to the 
arbitrariness of function u. One can assume the right-hand side of equation (18. 5p to be 
a definition of a new function u{k^, k^), i.e. n = 1, w = P = 0. At that, from the first 
equation of (18. 9p it follows that functions and are functionally dependent. The 
solution acquires the following canonical form: 



X = (A;i + F{T^{k\ e))) e' + {g{k')+T^{e, e))e^+T\k\ e)e' 



1.(8.11) 



Here F, and g are arbitrary functions. 

In case m 7^ by using the arbitrariness of function u{k^, k^), the right-hand side 
of the first of equations (18.50 can be taken as a new function u. Then m = 1, v = a = 0. 
Functions n, w and (3 in the second of equations (18.50 are changed correspondingly. 
Next, by using the equivalence transformations (12. 8p and (12.90 one makes w = 1, (3 = 0. 
Thus, 





= /(A)n(P) + 


k\ 


rl 


'2 


= -1, 




rl 


^3 







-n{k^ 



(8.12a) 



(8.125) 
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Finally, the solution takes the canonical form as 

X = (/(A) + T\k^, e))e^ + ((7(A) + T^{e, k^))e^ + T\k\ P)el (8.13) 

Here the relation f l8.12ap for arbitrary choice of functions /, g and n determines function 
\{k^^k^). Functions r^, r^, and should satisfy equations (18.126|) . 

Hence, equations (I7.2p are completely integrated. Solutions are given by formulae 
<Km and (EH. 

For s = 1 equations (17.31) are considered. Integration of the second of the equations 
(17. 3p leads to the following representation of vector cr: 

o- = M(A;\A;=^)e(P) + (5(P). (8.14) 

Here U = du/dk^. Differentiation of equation (I8.14p on variables fc^, k^ and substitution 
into the first of equations (17. 3p lead to the relation 

(MMie' + Mi5'- I/) X e = 0. (8.15) 

In order to separate variables in equation (I8.15p . cases of different space of the dimension 
of the linear space spanned by the 3D vector {uui,ui, 1) for various values of k^ should 
be investigated. By virtue of condition uu ^ the last two components of the vector 
are always linearly independent, hence the dimension can be either 2 or 3. 

i) dim{(uui, Ml, 1)} = 3. In this case all components are linearly independent, 
which allows the separation of variables in equation (I8.15P : 

e' X e = 0, <5' X e = 0, x e = 0. 

By the integration one finds 

e = e{k^)'n, 6 = S{k^)r], v = u{k^)r], ry = const. (8.16) 

Substitution of representations (I8.16p into the second series of equations (17. 3p leads to 
the contradiction. 

ii) dim{(uui, Ml, 1)} = 2. In this case there exists a linear dependence between 
components of vector (-u-ui, mi, 1). Due to the linear independence of components ui 
and 1 the linear relation can be chosen as 

uui + a{k^)ui + b{k^) = 0. 

By integration one finds 

^u^ + a{k^)u + b{k^)k^ = c{k^). 

Function u is determined accurate to an arbitrary additive that depends on k^. This 
arbitrariness allows one to zero function a. Next, the equivalence transformations (12. 8p 
and (12. 9p allows to normalize parameters as 6 = —1/2, c = 0. Thus, in this case function 
u is equivalent to the following one: 

u = ±y^\. (8.17) 
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Separation of variables in equation flS.lSp leads to the relation 

- e' - ly] X e = 0, 6' xe = 0. 
From this, accurate to insufficient additive, one obtains 

i/ = ^e', d' = a{k^)e. (8.18) 
Under this condition the second series of equations (17.31) reads 

e'-(eXT2) = 2, T3-(eXT2) = 0. (8.19) 

By virtue of the ffist of equations (I8.19P vectors e and T2 are linearly independent, 
whereas due to the second of this equations vector T3 can be linearly expressed in terms 
of T2 and e. This fixes vector r in the form 

T = ^{X)+u:{e), X = k^ + a{k^), Lj' = /3(P)e (8.20) 

with arbitrary vector $ and arbitrary functions a and f3. From the ffist of equations 
(I8.19P it follows the relation 

e'- (e X $'(A)) =2. (8.21) 
Independence of vector e on variable k'^ implies 

$ = Ar7 + ?7, C = const. 
Substitution into equation (I8.2ip yields 

e ■ {exT}) = 2. (8.22) 
Differentiation of equation (I8.22p on k^ gives 

e" . (e X 77) = 

this implies that vector e" can be expressed linearly in terms of e and t]. Hence, 

e = e^''''n^ +e-^''''n^ + 07) (8.23) 

with arbitrary vectors 7, C and arbitrary constant vectors rj^ and r}^. Substitution of 
representation fl8.23p into equation fl8.22p gives a restriction on vectors 77^ and 77^: 

777 ■ (t/^ X ry^) = 1. (8.24) 

Finally, accurate to insufficient terms, one obtains the solution 

X = ±^Ak^ (e^'^'ri^ + e'^^'rf + Cry) + ^'77+ 

r r r (8-25) 

+ / a{k^)e^'''dk^r]^ + / a{k^)e-^^'dk^r]'^ + C a{k^)dk^r]. 

Here 7 and C are arbitrary constant such that 7 may take complex values; in this case 
vectors r/* should be chosen such that the whole solution is real. Function a{k^) is 
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arbitrary; constant vectors 77, rj^, and rj^ satisfy equations fl8.24p . Thus, the canonical 
form of the solution of equations (17. 4p in the observed case is given by equation fl8.25p . 

Finally, for s = 1 equations fl7.4p are observed. By virtue of the second of equations 
fl7.4p vector cr depends linearly on and, hence, the solution is reduced to either of 
solutions (15. 2p or (15. 3p . This remark completes the investigation of case s = 1. 



9. Integration of equations in case s = 2 

Equations (17. 5p are now observed. By virtue of the arbitrariness in the choice of 
linearly independent functions and f/^ this system admits the following equivalence 
transformations: 

A^B''^A^B\ (9.1) 

K = K + uo^{k'^)e^ + uo'^{k^)e^. 

By the action of transformations (19.10 one can bring vectors e^, to the following form: 

el = el + a(P)e^ = + I3{e)e\ k = e^ (9.2) 

where a, /3 are arbitrary functions. 

Integration of the second equation of (17. 5p on k^ leads to the following representation 
for vector cr: 

cr = u\k\ e)e\e) + u\k\ k'^)e\e) + k^ii^e). (9.3) 

Here f/' = du^/dk^. An arbitrary additive vector depending on k^ that appears as a 
"constant of integration" can be zeroed by the suitable change of functions and v?, 
and by application of transformation (12. 9p . From representations (19. 2p . (19. Sp for vector 
cr one finds 

cr = u\k\ k^)e^ + u^{k\ k^)e^ + S{k\ k^)e\ S = au^ + (5u^ + k\ 

This implies 

ai = u\e^ +ule^ + Sie^, 0-3 = w^e^ + WgC^ + Ssel (9.4) 

Hence, 

(T^x ai = {u\u\ - u\ul)e^ - {u\Si - u\S:i)e^ + {ulSi - u\S:i)e^ = 
The latter equation is satisfied only if 

ulul = u\ul, ulSi = ulSs, ulSi = ulS^. (9.5) 
Equations (19.50 imply the functional dependence 

= u\X{k\k^)), = u\X{k\e)), S = S{X{k\e)). (9.6) 
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Function X{k^, k^) can be found from the following equation: 

au\X) + (3u\X) + k^ = S{X) (9.7) 

For some functions a{k^), (3{k^), u^{X), u^{X) and S{X). 

Finally, vector cr has the form defined by the equation (]9.3p with functions u^, 
u^, given by relations ( 19. 6p . ( 19. 7p and vectors e\ e^, k, satisfying equations ( 19. 2p with 
arbitrary functions a, p. 

Equations for vector T(fc^, k^) remain. After substitution of obtained representation 
for vector cr in equation (14. ip and splitting with respect to linearly independent functions 
and u"^ one obtains the following equations: 



(e^+ae^)-(T2XT3) 



0, 



(e2+/3e3).(r2XT3) = 0, e3.(T2XT3) = 1.(9.8) 



Vector T can be written in the basis of linearly independent vectors e^, e^, k as 
follows: 

T = ri(e^ + ae^) + r2(e2 + /3e^) + r^e^ (9.9) 

Here functions r* depends on fc^ and k^. Substitution of representation (19.91) in equation 
gives the following relations: 



Tn To 



0, 



(9.10) 



Tn To" 



(9.11) 



Tn To 



^3 ^3 



(9.12) 



By virtue of equation (I9.12p . equations (I9.10p and ( 19. lip have non-trivial solutions only 
if 



= 0, Tg^ + T^a' + T^l3' = 0. (9.13) 

As a consequence of equation ( I9.13p . function in representation (19. 9p depends only 
on k^. For simplicity the notation = 7(/c'^) is adopted. In the second equation of 
(I9.13P one distinguishes the following cases 

(i) a' = (3' = 0; 

(ii) {ar + + 0; 

In the first case 7' = 0, and functions r^, are such that equation (19.120 is satisfied. 
This leads to a solution 

X = {u^{k^) + ri(P, k^))e} + {u'ik^) + r2(P, P))e2 + fc^e^, (9.14) 
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where u , u are arbitrary functions; r , r are arbitrary functions satisfying equation 



^3 ^3 



1. 



In the second case by virtue of equation fl9.13p one obtains the representation 



a'y,{k\k')- — 



(9.15) 



Substitution of this representation into equation fl9.12p gives the following equation for 
function ip: 



V5T^(ttP -pa)- — la 



/3'7' 



1. 



9A;2^"" dk^ \ V (a')' + (/?') V V(a')' + 

Integration of this representation on k'^ and simplification of the second term leads to 
the quadratic equation for function ip: 



a'f3"-f3'a" , 
^ <^ 



7' 



„ 7' (a'a" + /3'/3") 



[a 



'\2 



^ = k\ 



(9.16) 



The final form of the solution is 

X = {u\\) + T^)e' + {u\\) + T^)e' + (^(A) + ar^ + I3t^ + 7)6^. (9.17) 

Here a{k^), f3{k^), lik'^) are arbitrary functions. Functions '^^(A), '^^(A), S{X) are also 
arbitrary and linearly independent; function \{k^, k^) is given by the implicit equation 
(19. 7p . Functions T'^{k'^, k'^) and r^(fc2, fc^) is specified by relation (]9.15p . in which function 
if is determined by the condition fl9.16p . 

For s = 2 in the case (17. 6p the system of equation admits the following equivalence 
transformations: 

= A\e)e^ + A^k'^^e^ 
_ A^B^^A^B\ (9.18) 

In the orthonormal basis of vectors e^, e^, e^, vectors e^, and r have the following 
components: 



.11 ^12 A'i\T 



.21 ^22 ^23\T 



-1 J2 ^3\T 



By the action of equivalence transformations (I9.18P vectors and can be brought to 
a form 

= e^ + a{e)e\ = + (3{e)e^ (9.19) 

where a and /3 are some arbitrary functions. 

First, equations for vector cr are observed. Substitution of relation fl9.19p into the 
second of equations (17. 6p gives 

<Ti = u\e^ + ule^ + {au\ + /3u?)e^ (9.20) 
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Integration of this equation on and differentiation on produces 

0-3 = ule^ + ule^ + {piu^ + au\ + /3V + + (9.21) 

Vector product of these formulas, collection of coefficients at basic vectors, and equating 
of corresponding coefficients in the right-hand side leads to the following system of 
equations: 

u\u\-u\ul = \, M}(a'M^ + /3V + y) = 0, M2(^aV + /3V + 7') = 0. (9.22) 

Two cases are possible. Case a) = = contradicts to the linear independence of 
functions = u\ and f/^ = m^. In case b) a'u^ + f3'u'^ + 7' = by virtue of linear 
independence of functions and U'^ one obtains a'{k^) = (3'{k^) = 'y'{k^) = 0. Then, 
functions a{k^), f3{k^) and 'j{k^) are constants, that can be taken zero without limiting 
the generality. 

Equations for r are now observed. Substitution of representations f l9.19p into 
equation for r gives the following relations 

7(fc3)r| = 1, (9.23) 



0, 



To To 



0. 



(9.24) 



Equations f l9.24p imply 

T = T{x{e,e)). 

By integration of equations fl9.23p on fc^ one obtains 

j(k^y = e 

The arbitrary constant of integration appearing at the integration can be zeroed by 
the equivalence transformation f l2.9p . By applying transformation (12. 8 p one can make 
7(fc^) = 1. By virtue of the equation for function one can assume A = fc^. 
The solution takes the following form: 

X = {u\k\ e) + T\e))e^ + {u^{k\ P) + T^{e))e^ + 



2„,1 



U^U 



1„,2 



1; function 



where functions u^{k^,k^), u'^{k^,k^) satisfy the relation 
r^(A;^), T'^{k'^) are arbitrary. This solution coincide with solution (I9.14p . 

In case of equations (17. 7p the following equivalence transformations are admitted: 

K = A{k^)K + B{k^)e, e = C{k^)e, A,C^O. 

Under this transformation vectors e, k, and r can be represented in the orthonormal 
basis of vectors e* as 



e = e' + a\e)e^ + a\e)e\ n = + I3e^ 



[T^ ,T ,^ ) ■ 



(9.25) 
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In the same manner equations for vector cr give 



u 



u\u^{a^)' + u\k^f5' -u\(5 = u^P, (9.26) 

ula" + {uW + l){u\a'y - k'(3') - - u\ay{u\a^ + /3) = {a^P - a^)u\ 

Two former equations of this system imply 

This equation is satisfied either if u\ = 0, or if u^{{a'^yf3 — (a^)') — k^f3' = 0. First case 
contradicts to the hnear independence of functions and U"^. The second case leads 
to equations 

- (a^)' = 0, /3' = 0. 

Hence, (3 = (3q = const, a^(3o — a^ = m = const. The third of equations (19.261) is satisfied 
identically. 

Now equations for vector r are observed. Substitution of representations (19.251) 
into equations (17. 7p gives 

rim-r|/3o + r| = 0, (9.27) 



Tn To 



TX To 



a 



To To 



Ti To 



a 



To To 



Ti To- 



(9.28) 



To To 



Tn To 



+ /3 



To To 



Ti Ti 



1. 



Form the first equation of the system it follows 

= S{k^) + (Sot"^ -mT\ 
Substitution of into equation (I9.29P yields 

T' = -J + U{k% 

By plugging functions r^, into equation (I9.28P one obtains 

2 2 

r = a 



The final form of the solution is the following: 



2\ / Pq;^\ 

■y]e^+ik^ + u^o? + -i + a^v — \ 



(9.29) 



(9.30) 



Here u^{k^,k^), n{k^), 6{k^), a'^{k^), 7(fc^) are arbitrary functions; m, Pq are arbitrary 
constants. 
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(9.31) 



Now equations fl7.8l) are investigated. This system admits the following equivalence 
t r ansf ormat ions : 

e = A{k^)e, V = B{k^)u + C{k^)^i + D{k^)e, 

Tl = E{e)^i + F{k^)e, A,B,E^O. 

Under the action of transformations fl9.3ip vectors e, fi and v can be represented in the 
orthonormal basis e* as follows 

e = + a\k^)e'^ + a^{k^)e^, = + f3{k^)e\ iy = e^. (9.32) 

Let us observe equations for cr. Substitution of representations fl9.32p into equation 
(17. 8p and comparison of respective components gives the following equations 

a^'i-ul + u\a^yul) + a\ul/3 - u\a^y + 1) = 0, 

-uj(3 + u\a^yu\ = 1, 

The first of equations is satisfied identically by virtue of the remaining two. Substitution 
of ul into the second equation of the system and integration with respect to k^ accurate 
to insufficient constant of integration gives an equation for function u^: 

{u'f{~(3ia^y + {a'y) = 2k\ 

Hence, 



1 



u 



2A;i 



(a3)' - /3(a2)' 



Vector cr takes the form 

(T = ± 



2k^ 



a^y-p{a^y 

Now we solve equations for vector r. Substitution of equations (I9.32p into (I7.8p 
gives the following system of equations: 

{a^ - a'l3)Tl + I3tI - r| = 0, -a\l + = 1, (9.33) 



^2 3 

^3 ^3 



— a 



^1 ^3 
^3 ^3 



a 



^3 ^3 



0. 



(9.34) 



Form the second equation of system (19.331) by integration on k"^ accurate to an insufficient 
function one obtains 



T = k + a T . 

Substitution of into the first equation of (I9.33P and integration on k'^ give 
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Substitution of and into the first equation of the system and integration yield 



k^(3' + y 



{a^y - /3(a2)' 
Finally, we obtain the solution in the form: 

2k^ k^P' + i 




- /3(a2)/ (a3)'_/5(a2)/ 





Here ^^(/c^), ^^(/c^), /^(/c^), 'l{k^) are arbitrary functions such that (a^)' 7^ 

It is remain to observe the last case given by equations f l7.9p . Similar calculations 
shows, that this system does not have non-trivial solutions. Thus, investigation of the 
case s = 2 is completed. 

10. The final list of solutions 

Analysis of systems f l7.ip - fl7.9p produces non-equivalent solutions given by the following 
formulae: ([EID, O, flHUj) . (l87[3|) . (ESS]), (EUD, dOrTTj) . (ICTj) and (lOSj) . At that, 
solutions f l8.25p . f l9.30p and fl9.35p represent special cases of formulae (15. 2p . (15. 3p . It is 
convenient to sort the obtained solutions in groups. 

According to the representation of solution (13. ip Maxwellian surfaces are 
translational, i.e. are sweeped by the parallel shift of the curve x = cr{k^,kQ) along 
another curve x = T(fc^, k^) provided k^ = k^ is fixed. Curves are parameterized by 
k^ and k"^. Each of these curves can be a straight line (zero curvature), planar (zero 
torsion), or a general 3D curve. For convenience the curve with the lower "dimension" 
will be called the generatrix, whereas the curve with greater "dimension" will be called 
the directrix. 

Theorem 2. Solutions of equation ^3.^) are exhausted by the following canonical 
representatives: 

(i) Solution with a straight generatrix, given in the general form by equations Ii5. 2\] . 

(El. 

(a) Solutions with a planar generatrix, given by formulae hS.ll]] . h8.1'J\] . ^9.1^^ . 
(Hi) Solutions with spatial generatrix, provided by formula ^9.11 ). 



The complete set of solutions is obtained from the canonical ones by application of 
equivalence transformations (dl^j and Ii2.9\) . 
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The following theorem introduce additional classification on the set of solutions. 

Theorem 3. All enumerated in Theorem IE solutions are either specializations of 
solution I{6.1a\) - f6.1c\) with appropriately chosen vectors a and cx, or solutions of the 
form Ii5. 3\) . 

Proof. For the proof it is required to demonstrate vectors cr and ex. in formulae fl6.1ap - 
(16 -Id . that provide the coincidences of solutions. 

1. Inclusion (15. 2p C (I6.1ap -( l6.1cp is provided by the following choice of vectors 

a- = {k\ 0, 0)^, cx = (1, a\e, k'), a\k\ k')f. 

Hence, the second and the third components of equation (I6.1cp can be taken as 
definitions of functions and a^. The first component of this equation gives the 
required restriction for functions and r^. 

2. Inclusion flSHD C fl833D . Let us take 

n(P) = 0, fiX) = k\ g{X) = g{k'). 
in equation (I8.13p . Accurate to the rotation of coordinate axis this gives he solution 

dHUD. 

3. Inclusion fl813ll C flO^ - flO^ . Let us take 

CT = {fiX),9iX), of, a = (-n(P), 1, aik', k')f. 

As before, an arbitrary function a{k'^,k^) is given by the last component of equation 
(I6.1cp . The remaining two components of the equation give the required restrictions to 
the vector r. 

4. Inclusion (191^ C flOoD - flOcD . Let us choose 

(T = {u\k^),u\k^),k'f, a = (0,0,1)^. 

Equation (16. Id) has a non-trivial solution only if = 0. The required equation for the 
remaining components of vector r follows from the third component of equation (I6.1cp . 

5. Inclusion figTTp C flO^ - flOcD . Here 

cr = {u\X),u\X), S{X)f, a = (-a(P), -P{k^), if. 

In this case equation (I6.1cp for vector r accurate to equivalence transformations has a 
solution, given by formulae (I9.13p . (I9.15p . (19.161) . Indeed, by scalar multiplication of 
equation (16. Id) on T2 and integration with respect to k"^ one obtains an integral 

r(fc2,A;2) ■a(A;3) = /(A;^). 

It implies, that vector r has the representation 

r = T\k\ k^)e\k'') + T\k\ k^)e\e) + H,{k^), 
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where 

e^-ct = 0, e^-a = 0, K-a = f{k^). (10.1) 

By virtue of arbitrariness in the choice of functions r^, r^, equivalence transformations 
(19. ip act on vectors e^, e^, and k,. Due to these transformations, taking into account 
equations (110. ip these vectors can be brought to a form (19. 2p . In this representation 
equation (I6.1cp is equivalent to the system of equations (19. 8p . All the remaining 
considerations coincide with those used above for obtaining of the solution. □ 

11. Conclusion 

It is shown that all solutions of ideal MHD equations, that describe stationary flows 
with constant total pressure belong to classes given by formulae (15. 3 p and (l6.1ap -( 1B.lcp . 
To prove this statement a curvilinear system of coordinates, in which streamlines and 
magnetic lines play a role of coordinate curves was used. It was shown, that Maxwellian 
surfaces, spanned by streamlines and magnetic lines, belong to a class of translational 
surfaces. Formulae (15. 3p and (l6.1ap -( 1B.lcp provide restrictions of possible forms of 
directrix and generatrix of the surfaces. 

Maxwellian surfaces determined by solution (15. 3p are cylindrical. The directrix can 
be chosen with functional arbitrariness according to formulae (15. 3p . (15. 4p . The direction 
of the generatrix can change while moving from one Maxwellian surface to another 
according to the modification of vector r]{k^). Functional arbitrariness allows one to 
control the shape of streamlines and magnetic lines on Maxwellian surfaces. 

For the family of solutions described by formulae (I6.1ap - (l6.1cp all Maxwellian 
surfaces have common directrix, determined by the dependence cr(A). Functional 
arbitrariness of the solution allows varying of the picture of the described fluid motion. 
In particular, solutions of this type describe stationary flows in arched jets, that can be 
used for description of solar prominences. For practical computations it is convenient 
to use special forms of solutions determined by formulae (18. lip . (I9.14p . 

The work is supported by RFBR (grant 11-01-00026), Presidential programs N.Sh- 
4368.2010.1 and MD-168.2011.1, and Federal Purpose Program 02.740.11.0617. 

Appendix A. Separation of variables in equation ( 13. 2 j) 

According to Lemma [1] and by virtue of symmetry of equation it is sufficient to observe 
cases s = 0, 1, 2. 

Case g = 0. Here matrix C is absent. Let us express the 6-dimensional vector d as 
a union of two 3-dimensional: d = {fi{k^), u{k^)). Hence, 

0-3 X O-i = ^(P), (Ti=l^{k^). (A.l) 
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At that, in order to guarantee the dependence of vector cr on coordinate k^, it is supposed 
that u ^ 0. Note that accurate to equivalence transformation, this solution belong to 
class of solutions with cylindrical Maxwellian surfaces, described in section |5l However, 
for the sake of completeness, this case will be observed separately. 

The scalar multiplication of left- and right-hand sides of equations f lA.ip gives 
fi ■ u = 0. This equation allows representation of vector fj, 

fl = fj, X u 

with some vector fi{k^). By plugging the obtained representation of vector into the 
first equation of (IA.1|) by virtue of the second equation one obtains 

(T:i X ai = ft X u = ft X (Ti =^ (erg — p.{k^)) x cTi = 0. 

Accurate to equivalence transformations one can take fj-ik^) = 0. Thus, cr^ x cri = 0, 
and system of equations (17.11) follows. 

Case g = 1. Lemma [1] gives the following equations 

(T3X (Ti = U{k\k^)^i{k^) + u{k^), (Ti = U{k\k^)e{k^) + K{k^) (A.2) 

with some auxiliary function U {k^, k^) such that dU/ dk^ 7^ 0, and 6-dimensional linearly 
independent vectors (/x, e) and {iy,K,). Cases of different dimension of the linear space 
{e, K,} spanned by vectors e and k for various values of parameter k^ are observed. 

1) dim{e, k} = 2. By scalar multiplication of left- and right-hand sides of equations 
f lA.2p and by splitting on the different powers of function u one obtains the following 
relations: 

fi ■ e = u ■ K, = 0, fi ■ K, -\- 1/ ■ e = 0. (A. 3) 

From the former two equations of system (1A.3P the representations follow: 

^ = fj, X e, v = i> X K (A. 4) 

with suitable vectors fi and u. By plugging this equalities into the latter equation of 
(lA.Sp one obtains 

{fi X e) ■ n + {u X k) ■ e = Q =^ (i> — /i) ■ (e x k) = 0. 

Hence, 

u = fL + a{k^)e + l3{k^)K (A.5) 

with some functions a, /3. Since vectors ft and u are defined accurate to additives, 
proportional to vectors e and k, respectively, in equations (lA.Sp without limiting of 
generality one can assume a = P = 0. By using representations f lA.4p with u = jl in the 
first equation of (lA.2p . taking into account the second equation of (lA.2p . and collecting 
terms one obtains 

(0-3 - fi) X ai = 0. 
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By virtue of dependence of vector fi only on and due to the arbitrariness of vector cr, 
without hmiting the generahty, one can assume fi = 0, i.e. 0-3 x cti = 0. Substitution 
of obtained representations into equations (14. ip and sphtting on powers of function U 
gives equations for vector r, specified in system (17. 2p . 

2) dim{e, At} = 1. Two cases are possible: either a) k = a{k^)e, or b) e = and 
K ^ 0. In case a) by a suitable choice of function U in representation (1A.2P one can 
make k = 0, i.e. 

(Ti = U{k\k^)e. (A.6) 

Dot product of obtained equation with the first of equations (IA.2p . and splitting with 
respect to powers of function U we obtain equalities, analogous to equations (1A.3P : 

^■€ = 0, iy-e = 0. (A.7) 

Equations (1A.7P yield the representation 

= fi X e, V = i> X e (A. 8) 

with suitable vectors fi and u. Substitution of representations (lA.SP into the first of 
equations flA.2p taking into account equation flA.6P leads to the relation {0-3 — fi) x cti = 
z> X e. Arbitrariness in the choice of vector cr allows simplification fi = 0. Thus, for 
vector cr one obtains the following relations: 

erg X cri = j> X e, (Ti = U{k^ , k^)e. (A. 9) 

Substitution of equation flA.9p into (14.11) and splitting with respect to U results equations 
(17. 3p for vectors cr and r. 

Similar actions in case b) leads to equations (17.41) . 

Case s = 2. The following equations are to be observed: 

0-3 X cri = U\k\ e)fi\k'') + U^{k\ k^)n\k^) + iy{e), (A.IO) 

= U\k\ k'')e\e) + U\k\ e)e^{e) + k(A;3) (A.ll) 

Here 6-dimensional vectors (^^,e^), (/x^,e^), (i^, k) are linearly independent for every 
value of parameter k^; functions f/^, U"^ and 1 are linearly independent as functions of 
k^ (coefficients in linear combinations of functions can depend on k^). The classifying 
parameter is the dimension of the linear space {e^, e^, k} spanned by vectors e^, e^, and 
K for different values of parameter k^. 

1) dim{e^,e^,K} = 3. The dot product of left and right-hand sides of equations 
flA.lOp . and flA.lip . sphtting with respect to linearly independent functions f/^, f/^, and 
1 yield: 

/x^ ■ = • = 1/ ■ K = 0, 

(A.12) 

/x^ ■ + ^2 • = 0, • K + J/ ■ = 0, /x^ ■ k + • = 0. 
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From the former three equations of (1A.12P the representations follow: 

fj,^ = fj,^ X e\ /x^ = /i^ X e^, u = x k (A.13) 

with some vectors u. Substitution of these representations into the latter three 
equations of (1A.12P gives 

(A^ - A^) ■ (e^ X e^) = 0, (A^ - u) ■ (e^ x k) = 0, (^^ _ . ^ k) = 

Taking into account linear independence of vectors e^, e^, k, and the arbitrariness of 
vectors jl^, u, this implies 

fj,^ = fi^ = V. 

Substitution of formulae (lA.lSp with the preceding equations into ( 1A.10I) after some 
simplification leads to 

0-3 X CTi = i> X {U^e^ + f/^e^ + k) = x ai =^ {a^ — i>) x cri = 0. 

Since vector cr is defined accurate to the additive vector, depending only on k^, without 
restriction of generality one can assume i> = 0. 

Substitution of the obtained representation of vector cr into equation (14. ip and 
splitting with respect to the remaining linearly independent functions U^, U"^, and 1 
produces equations (17. 5p . 

2) dim{e^,e^,K} = 2. The following subcases are possible: a) vector k, can be 
linearly expressed in terms of vectors and e^, or b) vectors and are proportional 
or one of them is equal to zero. 

In case a) by addition of terms, depending only on to and U"^, one can make 
K, = 0. Dot product of left- and right-hand sides of equations flA.lOp . flA.lip and splitting 
with respect to functions f/* leads to equations (1A.12P with k = 0. These equations, in 
turn, give 

= fix e^, /^^ = /X X 1/ = 6{k^)e^ x (A.14) 

with some vector fiijc") and function dijc"). Substitution of obtained representations 
into equation flA.lOp and action of equivalence transformations lead to equations (17.60 . 

In case b) by a proper non-degenerate change of functions [/^, and f/^ it is always 
possible to make = in equation flA.lip . In this case relations (IA.12p produce 
representations for vectors and u in the form 

fj,^ = fx X e^, fi^ = ae^ x k, u = fj, x k (A. 15) 

with some vector fi. The system of equations (17. 7p is obtained in the same way as in 
previous case. 

3) dim{e^,e^,K} = 1. The following subcases are possible: a) vectors k and 
linearly related with vector e^, or b) vectors are zero: = = 0. By the analogous 
calculations one can show that case a) gives a system of equations (17. 8p . and case b) 
reduces to the system (17. 9p . 
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